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ABSTRACT 

We consider the possibility that massive primordial black holes are the dominant form of dark matter. 
Black hole formation generates entropy fluctuations that adds a Poisson noise to the matter power 
spectrum. We use Ly-a forest observations to constrain this Poisson term in matter power spectrum, 
then we constrain the mass of black holes to be less than ~ few x 10 4 M Q . We also find that structures 
with less than ~ 10 3 primordial black holes evaporate by now. 

Subject headings: Cosmology: theory — dark matter — black hole physics — large-scale structure of 
universe — cosmology: observations — quasars: absorption lines 



1. INTRODUCTION 

The nature of Cold Dark Matter (CDM) is one of the 
most important and profound mysteries in modern cos- 
mology. One of the clues towards solving this mystery is 
the fact that Big Bag Nucleosynthesis (BBN) paradigm 
(e.g., Buries et al. 2001) predicts a much lower density 
for the baryonic matter (by a factor of 5 — 10) than al- 
most all the measurements of the CDM density in the last 
decade (see Bahcall et al. 1999 for a review). Black holes 
are an alternative to the most popular candidate for the 
dark matter: massive, non-baryonic, elementary particles 
(see Khalil & Munoz 2002 for a recent review). If Pri- 
mordial Black Holcs(PBHs), formed in the early universe, 
they could become the dominant form of dark matter. 

There are variety of mechanisms proposed for PBH for- 
mation: during inflationary reheating (Garcia-Bellido & 
Lindc 1998), from a blue spectrum of primordial density 
fluctuations (Carr & Lidsey 1993) or during a phase tran- 
sition in the early universe (e.g., Rubin et al. 2000). The 
typical mass of PBHs at the formation time can be as 
large the mass contained in the Hubble volume, Mjj, down 
to around 10" 4 Mh (Hawke & Stewart 2002). Carr & 
Hawking (1974) show that there is no physical self-similar 
growth mode for the PBHs in a radiation gas, and so there 
will be no significant accretion after PBH formation. How- 
ever, Bean & Magueijo (2002) speculate that the accretion 
of quintessence (i.e., a scalar field) into the PBHs could 
lead to the scaling of the PBH mass with Mr. 

In this letter, we consider PBHs with 10 M < M PBH < 
10 6 M as the constituents of CDM. The dynamical con- 
straints on such objects has been thoroughly reviewed and 
discussed by Carr & Sakellariadou (1999). They show that 
the strongest upper limit on the mass of such compact ob- 
jects in the galactic halo is ~ 10 4 M Q which comes from 
two independent methods: (1) the rate of globular cluster 
disruptions and, (2) the mass of the galactic nucleus, ac- 
cumulated due to dynamical friction. However, the glob- 
ular cluster disruption limit depends on the completely 
unknown initial number of globular clusters, while the 
galactic nucleus limit ignores black hole ejection through 
gravitational slingshot which may prevent mass accumu- 



lation at the center (Xu & Ostriker 1994). Therefore we 
will take a more conservative value of 10 6 Mq as our up- 
per limit. On the lower end, the non-detection of long 
duration lensing events by the MACHO microlensing ex- 
periment excludes the mass range 0.3 — 30 M Q for the halo 
objects (Alcock et al. 2001). This letter studies the impact 
of these PBHs on the large scale structure of the universe, 
and specifically the signatures of discreteness in CDM. In 
Sec. 2 we study the impact of the Poisson noise due to the 
discreteness of PBHs on the linear power spectrum. Sec. 3 
compares the numerical simulations with this linear power 
spectrum with the observational data of the Ly-a forest 
power spectrum. This puts an upper limit on the mass 
of PBHs. In Sec. 4 we discuss the effect of evaporation 
of early structures, which puts a lower cutoff on the mass 
function of the PBH haloes. Finally, Sec. 5 concludes the 
paper. 

2. THE LINEAR SPECTRUM 

In the following analysis we assume that PBHs form as 
the result of a phase transition at some temperature T c , 
in the radiation dominated era. Furthermore, we will as- 
sume that there is no substantial accretion following their 
formation era and all of their masses are around ~ M PBH . 
This assumption can be generalized to an extended accre- 
tion scenario where T c marks the end of accretion period 
(Magueijo and Bean 2002). The black hole density follows 
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Sp is the fluc- 



where no, PBH is the number density of PBHs, 
tuation is the number of PBHs due to Poisson noise, and 
the T/T c factor describes the dilution in PBH density due 
to cosmic expansion. Since the PBHs are not correlated on 
acausal distances, one expects that on scales larger than 
the Hubble radius 
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For linear perturbations, on large scales, Eq.(l) leads to 
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where the 5 rBH ,S p and S r are the relative overdensities 
of PBHs, Poisson fluctuations and radiation, respectively. 
Since S p in Eq.(l)is observable and constant, one would 
conclude that the quantity 

3 

S = S pbh ~ = S p (4) 

is gauge-invariant and conserved. Indeed this is the en- 
tropy per PBH, which should remain constant as long as 
the universe expands adiabatically (e.g. see Mukhanov 
et al. 1992). The associated perturbations, generated in 
this way are isocurvature(or entropy) perturbations, as the 
curvature at large scales is not (immediately) affected by 
the formation of compact objects at small scale. 

As we are assuming that PBHs are the present day Cold 
Dark Matter (CDM), the overdensity of CDM is given by 

i CBM (fc) - T ad (k)5 i>ad (k) + T iso {k)S(k), (5) 

where T adi (k) and Ti SO (k) are the transfer functions for 
adiabatic and isocurvature perturbations respectively. For 
the following analysis we will use the analytical fits quoted 
in Bardeen etal. 1986 to the transfer functions. Eq. (5) 
leads to the following power spectrum 

P CD ,M = T 2 ad (k)P^ ad {k) + T? so (k)P p . (6) 

In this expression, P^ ad {k) — Ak n with n ~ 1 is the adia- 
batic power spectrum which is produced through inflation 
(or an alternative method of generating scale- invariant adi- 
abatic perturbations), while P p is given in Eq.(2). 

One can easily see that the isocurvature term on the 
RHS of Eq.(2) contributes a constant to the power spec- 
trum as both P v and 



T iso {k) = -(1 + z eq ) forfc > a eq H eq (7) 

are independent of k (e.g. Peacock 1998). Note that this 
is the simple linear growth due to gravitational cluster- 
ing which is the same for adiabatic fluctuation. Since the 
power spectrum of adiabatic fluctuations decays as k~ 3 at 
small scales, one expects to see the signature of this Pois- 
son noise at large fc's. Combining Eqs. (2), (6) and (7) 
gives the power offset 

Ap „ 9Af PBH (l + z eq f 
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Fig. 1. — Linear power spectrum for different masses of the PBHs. 
cTg is us f° r the model without the PBHs and the amplitude of the 
(initially) adiabatic modes is the same for all models. 

Fig.(l) shows the linear power spectrum for different 

masses of the PBHs. We see the Poisson plateau (Eq. 
8) at large k's which drops with decreasing mass. The 
impact of this plateau on the Ly-a forest power spectrum 
is discussed in the next section. 
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Fig. 2. — Influence of PBHs on the Ly-a forest flux power spec- 
trum, Pp(k). The black, solid curve shows our prediction for Pp(k) 
in a standard ACDM model (i.e., no PBHs) in which the amplitude 
of the linear power spectrum, <Tg, was adjusted to match the data 
points from Croft et al. (2002). The other curves show the predicted 
Pi? (A:) when white noise power due to PBHs with various masses is 
added. The Ly-a forest model parameters and Og were not adjusted 
to find a best fit for each mass so the disagreement between the PBH 
models and the data points does not indicate that the models are 
ruled out. 



which is also a lower bound on the matter linear power 
spectrum. 



3. SIMULATIONS OF LY-a FOREST 
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The lines in Fig. (2) show the predicted change in the 
power spectrum of the Ly-a forest transmitted flux, Pp(fc), 
as M PBH is varied. The points with error bars are Pf(&) 
measured by Croft et al. (2002) using their fiducial sam- 
ple (z — 2.72). The predictions were made using the large 
set of numerical simulations and the interpolation code de- 
scribed in McDonald et al. (2003). We have not yet per- 
formed fully hydrodynamic simulations using PBH power 
spectra, so our result is based entirely on hydro-PM sim- 
ulations (e.g., Gnedin & Hui 1998; McDonald, Miralda- 
Escude, & Cen 2002; McDonald 2003). The curves we 
show are smooth because the power spectra computed 
from the simulations have been compressed into the pa- 
rameters of an analytic fitting formula. The background 
cosmological model used in Fig. (2) is assumed to be flat 
with a cosmological constant, £1 CDM — 0.26, Sib = 0.04, 
h = 0.72, and n — 0.9 (this value of n is close to the best 
fit found by Croft et al. 2002 for our model). The Ly-a 
forest model assumed in the simulation is controlled by 
3 parameters: the mean transmitted flux fraction in the 
forest, F, and the parameters, T1.4 and 7 — 1, of a power- 
law temperature-density relation for the gas in the IGM, 
T = Ti. 4 (A/1.4)' 7 - 1 , where A is the density of the gas in 
units of the mean density (see McDonald 2003 for a demon- 
stration of the effects of these parameters on the flux power 
spectrum). The allowed range of each of these parameters 
has been constrained by independent observations. We 
use the measurement F — 0.746 ± 0.018 from McDonald 
et al. (2000) and the measurements T 1A = 20500 ± 2600 
K and 7 - 1 = 0.4 ± 0.2 from McDonald et al. (2001). To 
obtain these values at z = 2.72, we interpolated between 
the redshift bins used by McDonald et al. (2000, 2001). 
We subtracted 50% of the potential continuum fitting bias 
they discuss from F, and add the same number in quadra- 
ture to their error on F. We add 2000 K in quadrature to 
the error bars on T1.4 to help absorb any systematic errors. 
To produce Fig. (2), we fixed these Ly-a forest parame- 
ters to their measured values, and fixed the normalization 
of the initially adiabatic component of the linear power 
spectrum, <7g, to the value that gives the best fit when 
M PBH = 0. It is not surprising to see that the Ly-a for- 
est power increases dramatically as the white noise power 
from the PBHs becomes significant on the observed scales. 

Fig. (2) is not sufficient to place constraints on M PBH 
because we have not varied any of the other parameters to 
see if the predicted power can be adjusted to match the 
observations. To obtain an upper limit on the PBH mass 
we compute X 2 (Af PBH ), minimizing over the amplitude of 
the linear power and the three Ly-a forest parameters, 
subject to the observational constraints described above 
on F, T1.4, and 7 - 1. We follow Croft et al. (2002) 
in using only P^(fc) points with k < 0.04 s/km. Defin- 
ing an upper limit by X 2 (M r - Bli ) — X 2 (0) = 4, we find 
M PBH < 17000 M s . Fig. (3) shows how this limit is ob- 
tained. The temperature-density relation parameters play 
no role, but F is important. As M PBH increases, the best 
fit value of F also increases until this trend is halted by 
the external constraint on F. The effect of the increase in 
F is to reduce Pp(fc) (McDonald 2003), counteracting the 
increase in power due to M PBH . Finally, the unconstrained 
parameter Cg also increases with M PBH , further decreas- 
ing the power on small scales while increasing Pf(^) on 



large scales (see McDonald 2003; this freedom to adjust 
the tilt of Pp(fc) is significant - with Og fixed we find 
M PBH < 5800 M Q ). 
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Fig. 3.— Joint fits to the observed Ly-« forest mean absorp- 
tion level, F, and power spectrum, Pp(k), for several values of the 
PBH mass. Curves show the best predicted Pp(fe) for each value 



-'PBH ' 

F, T1.4, and 7 



1. The best fit values of F for the different 



values of M PBH are shown by the labeled horizontal lines. The large 
square is the measured value of F from McDonald et al. (2000). 

Fig. (3) may at first appear unconvincing to the reader 
unfamiliar with the Ly-a forest; however, the result is ul- 
timately simple to understand. In Fig. (1) we see that 
the white noise power begins to dominate on the scales to 
which the Ly-a forest is sensitive when M PBH ~ 10000M Q 
(note that 1 comoving Mpc/h =108 km/s at z — 2.72 in 
our model). As M PBH increases there is simply too much 
power on the scale of the Ly-a forest to produce the ob- 
served level of fluctuations. Increasing F can cancel some 
of the effect but the size of the increase is limited because 
F is directly observable. 

A factor of ~ 2 relaxation in the upper bound seems un- 
likely but not inconceivable. For example, if we arbitrarily 
increase the error bar on F to ±0.03, the limit we derive 
is M PBH < 41OOOM . The limit is M PBH < 37OOOAf if 
we arbitrarily decrease the predicted Pp(fc) by 10% for all 
models. The assumed value of n has no effect on the re- 
sult (we obtain M PBH < 18000Af Q using n — 1). Finally, 
we remind the reader that the Ly-a forest only constrains 
the power spectrum in km/s units at z — 2.72. Equa- 
tion (8) and our assumed cosmological model were used to 
compute M PBH . 

4. EARLY STRUCTURES AND RELAXATION EFFECTS 

So, apart from the Poisson noise, is there any difference 
between the gravitational clustering of the conventional 
CDM (WIMP particles) and PBHs? The answer is yes. 
The collisional relaxation time for a gravitational system 
is of the order of the number of particles, times the dy- 
namical time of the system. Therefore, one expects the 
relaxation related effects, e.g. evaporation and core col- 
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lapse, to happen faster for systems with smaller number 
of particles. As the structures form bottom-up in a hier- 
archical structure formation scenario (and even more so in 
the presence of PBHs as the spectrum is bluer), and the 
dynamical time for cosmological haloes is of the order of 
the age of the universe, such effects may be important only 
for the first structures, that form right after the recombi- 
nation era, which have the smallest number of PBHs. 

Let us make a simple estimate of how evaporation of 
early structures sets a lower limit on the mass of virializcd 
objects. The evaporation time of an isolated cluster can 
be estimated using 
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(see Binney & Tremaine 1987), where the subscripts, 
"evap" and "rel" refer to the characteristic times, sub- 
sequently associated with the evaporation and relaxation 
of the structure, while N and are the number of parti- 
cles and the median radius, respectively. To relate to 
the formation time, for simplicity, we consider a truncated 
singular isothermal sphere which within the spherical col- 
lapse model (e.g., Afshordi & Cen 2002) yields 
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where tf is the formation time of the object. Combining 
this with Eq.(9) gives 
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The next assumption is that the approximate formation 
time of the structure is when the variance of fluctuations 
at the mass scale of the structure a(M) is around the crit- 
ical overdensity in the spherical collapse model, 5 C (Gunn 
& Gott 1972), 

^W| form ~5 C ~1.67. (12) 

Based on the calculations of Sec. 2, a(M), which is dom- 
inated by the Poisson noise at the small scales, is given by 
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neglecting a late time cosmological constant which is not 
important at the formation time of the early structures. 
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Now, combining Eqs. (11-13) with 

1 ~ n o.2 ( !+£/ 
tf~ m \ l + z 

for a flat universe, gives the minimum mass for the struc- 
ture not to evaporate 
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* 3 x 10 3 M PBH (1 + z)- 1 (^) (A) 12/7 , (15) 

and consequently, the structures with M < M m [ n should 
evaporate by redshift z. 

This calculation may be an underestimate of M m i n since 
Eq.(9) is valid for an isolated object and ignores tidal ef- 
fects. In general, tidal stripping can decrease t evap signif- 
icantly (by a factor of 5-10). However, accretion, on the 
other hand, can slow down the net mass loss rate of the 
structure and hence has a reverse effect. In reality, the ex- 
act lower limit will probably depend on the environment 
of the structures. 



5. CONCLUSIONS 

We have studied the possibility of having CDM in the 
form of PBHs (or any other massive compact object) and 
its impact on the large scale structure of the universe. We 
see that the simple Poisson noise, enhanced by gravita- 
tional clustering in the matter dominated era, leads to a 
plateau in the power spectrum at large wave numbers (see 
Fig. 1). Comparison of numerical simulations of Ly-a for- 
est with the current observational data rules out the PBH 
masses larger than few x 1O 4 M . The discrete nature of 
the PBHs can also lead to the evaporation of small (early) 
structures. A simple estimate shows that this puts a lower 
limit of about 10 3 M PBH on the mass of small structures. 

PM thanks Nick Gnedin for the hydro-PM code. Also, 
NA and DNS are grateful to J.P.Ostriker for critical and 
illuminating discussions. The Ly-a forest simulations were 
performed at the National Center for Supercomputing Ap- 
plications. 
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